Freezing by Heating in A Driven Mesoscopic System 
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We investigate a simple model corresponding to particles driven in opposite directions and interacting 
via a repulsive potential. The particles move off-lattice on a periodic strip and are subject to random 
forces as well. We show that this model — which can be considered as a continuum version of some 
driven diffusive systems — exhibits a paradoxial, new kind of transition called here "freezing by 
heating". One interesting feature of this transition is that a crystallized state with a higher total 
energy is obtained from a fluid state by increasing the amount of fluctuations. 



Most of the phenomena in our natural environment 
occur under far-from equilibrium conditions resulting in a 
rich behavior both in time and space. An important class 
of such processes takes place in so-called driven systems, 
which have attracted considerable interest recently. In 
many of these systems, particles are driven either by an 
external field (force) , or they are self-propelled j|,D , 
and their collective behavior manifests itself in new kinds 
of transitions, including noise induced ordering ||J|] or an 
ordering in a continuous 2d velocity space 

Phase transitions are also common in equilibrium sys- 
tems, and the related analogies have represented an 
important contribution to the understanding of non- 
equilibrium processes. In most cases, lattice models have 
been considered to demonstrate non-equilibrium transi- 
tions. For example, jamming transitions have been seen 
in discretized traffic models |tJ and driven lattice gases 
P|. However, off-lattice (i.e., continuum) symmetry is 
known to bring in qualitatively new behaviour; in partic- 
ular, this is definitely so in 2d, see, e.g., the XY model 
versus the Ising model (in equilibrium). A continuum 
model may lead to new effects due to the fact that the 
notions of order and disorder have extra facets in this 
case, and it can describe compressible systems in a more 
delicate way. 

In this paper we will consider a simple continuum 
model exhibiting a paradoxial, new kind of transition 
that we call "freezing by heating" and being closely re- 
lated to situations relevant from the practical point of 
view. The model consists of particles driven in oppo- 
site directions and interacting through a simple repulsive 
potential. The particles move off-lattice on a periodic 
strip (in a two-dimensional tube) and are subject to ran- 
dom forces as well. The most interesting feature of the 
transition we find for this system is, that a crystallized 
state with a higher total energy is achieved from a fluid 
state over a transient disordered state by increasing the 
amount of fluctuations. 

In addition to the interest in the properties of driven 
systems on its own, there are several further motivations 



to study such a model. A system of light (rising) and 
heavy (sinking) particles in a vertical column of fluid, 
pedestrians moving in a passage, or a system of oppo- 
sitely charged colloidal particles in an electric field repre- 
sent potential applications of our model. In fact, the sys- 
tem we study is a generalization to the continuum case of 
a two-species driven lattice gas model proposed recently 
]l(| , with a number of relevant modifications arising from 
the adaptation to the off-lattice case. In a wider context, 
these models can be considered as simplified paradigms 
of systems consisting of entities with opposing interests 
(drives) . In the present work, we consider the behavior of 
a limited number of particles in a confined geometry, and 
our results are primarily valid for this "mesoscopic" sit- 
uation. In the quickly growing literature on mesoscopic 
systems there are many examples of the potential practi- 
cal relevance of phenomena occurring in various models 
for finite sizes |Tl| ]. 

We denote the location of particle i at time t by Xi(t) 
and its velocity dxi(t)/dt by Uj(t). Furthermore, we as- 
sume the acceleration equation 

dvi(t) voe, - Vi(t) 

m— — = m b&(f) 

at t 

+ + /b (*<(*)) . (i) 

to is the mass of the particle, vq the velocity with which 
it tends to move in the absence of interactions, r the cor- 
responding relaxation time, and £ {(1, 0), (—1, 0)} the 
direction into which particle i is driven. 7 = m/r may 
be interpreted as a friction coefficient, fij represents the 
repulsive interactions between particles i and j, /b the 
interactions with the boundaries, and £j the fluctuations 
of the individual velocities. For the interactions between 
the particles, we have chosen the simple function 

f ij (x i ,x j ) = -VA(d ij -D)- B , (2) 

depending on the parameters A and B, and the distance 
d%j[t) — \\ x i(t) ~ x j{t)\\ > D only. Thus, fij describes 
the effect of a soft repulsive potential of particle j with 
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a hard core of diameter D, reflecting the space occupied 
by the particle. Our choice of these details of the model 
corresponds to a motion of finite sized particles tending 
to avoid collisions and maintaining, if possible, a given 
velocity vq. In addition, the interactions with the bound- 
aries were assumed to be 

f h ( Xi ) = -VA(d iX - D/2y B , (3) 

where di± denotes the shortest distance to the closest 
wall. 

In contrast to previous studies of similar models Q , 
we will now investigate the decisive role of the fluctua- 
tions £i(t), which have been assumed to be uncorrelated 
and distributed according to a truncated normal distribu- 
tion with vanishing mean value and finite variance 9 [fl3| . 
We started our simulations with N particles randomly 
distributed on a strip without allowing overlaps. For half 
of the particles a driving into the (—1,0) direction, and 
for the other half a driving into the (1,0) direction was 
assigned. Numerical integration of equation (|l|), using 
periodic boundary conditions, has produced the follow- 
ing results: For small noise amplitudes 9 and sufficiently 
small particle densities, our simulations lead, depending 
on the strip width and the initial condition, to the forma- 
tion of two or more coherently moving linear structures 
(just as if the particles moved along traffic lanes) (cf. 
Fig. |l|a) . For relatively large N (if the available area is 
too small to allow freely moving lanes), jamming occurs. 
For a small intermediate density region, we find lane for- 
mation or jamming, depending on the respective initial 
condition. 

At small noise amplitudes, the mechanism of lane for- 
mation, which produces a "fluid" state, is very dominant 
and robust in our model. This can be understood as fol- 
lows: Particles moving against the stream or in areas of 
mixed directions of motion will have frequent and strong 
interactions, because of high relative velocities. In each 
interaction, the encountering particles move a little aside 
to pass each other. This sidewards movement tends to 
separate oppositely moving particles. Nevertheless, jam- 
ming may sometimes occur, but in most cases it also 
supports lane formation (see next paragraph). Particles 
moving in uniform lanes have very rare and weak interac- 
tions. Hence, the tendency to break up existing lanes is 
negligible, when the fluctuations are small. Furthermore, 
the most stable configuration corresponds to a state with 
a minimal interaction rate and is related with a maximum 
efficiency of motion [Q . 

Whereas spontaneous lane formation was also observed 
in previous studies of related models with deterministic 
dynamics only |lSq j, in the present, more realistic model 
we have discovered a surprising phenomenon when we in- 
creased the noise amplitude. If the fluctuations and the 
particle number are large enough, the particles crystal- 
lize into a hexagonal lattice. This is a consequence of 



several subsequent steps: First, the fluctuations are able 
to prevent lane formation or even to destroy previously 
existing lanes. This is so, because sufficently strong dif- 
fusion can prevent structure formation. Second, some 
of the oppositely moving particles block each other lo- 
cally from time to time. Third, this gives rise to jam- 
ming since, meanwhile, additional particles arrive at the 
boundaries of the blocked area. Fourth, if the jam exists 
long enough, both of its ends expand over the full width 
of the strip and develop "flat" boundaries perpendicu- 
lar to ei, in order to reach a balance of forces. For the 
same reason, the particles tend to arrange in a hexagonal 
lattice structure, very much like in a crystal. Fifth, the 
crystal is only stationary, if also the interface between the 
oppositely moving particles is, by chance, flat enough (cf. 
Fig. |l|c). In most cases, however, the interface is rough 
(cf. Fig. |l]b), i.e. in some of the horizontal layers, a 
majority of particles is pushing in one direction. As a 
consequence, the most advanced part(s) of the interface 
eventually break(s) through, which requires a continuous 
model, where the distance kept among the particles is 
flexible enough. In this way, particles with uniform direc- 
tions of motion form "channels" , which tend to produce 
lanes at sufficiently small densities and noise intensities, 
otherwise the particles jam again and again (as described 
above), until they end up in a stationary crystal. 

Due to the above described mechanism the crystallized 
state is metastable, i.e., sensitive to structural perturba- 
tions (like the interchange of a few particles in our case) . 
The crystallized state can also be destroyed by ongoing 
fluctuations with extreme noise amplitudes giving rise to 
a third, disordered ("gaseous") state with randomly dis- 
tributed particles. Thus, with increasing "temperature" 
9, we have the untypical sequence of transitions fluid — > 
solid — * gaseous. 

Interestingly, for a range of moderate densities we find 
a "fluid" state with lanes at small noise amplitudes most 
of the time, but a crystallized ("frozen" ) state, if the noise 
amplitude ("temperature") is large. We call this transi- 
tion "freezing by heating" . Starting with random initial 
conditions, the transition is rather smooth (cf. Fig. ||). 
This is partially so because the system can also become 
frozen at a relatively low noise amplitude, if the disorder 
in the initial state (in the sense of the deviation from a 
freely moving lane state) is large enough (which has an 
effect similar to additional fluctuations). The transition 
becomes sharper, if we always start with a two-lane state 
but with different random seeds. In any case, the transi- 
tion is hysteretic, since the noise-induced "frozen" state 
remains, when the noise amplitude is reduced, again. 

To characterize the state of the system, we calculated 
various quantities. The expression 

E= lim A Lif^, (4) 

l ~ l 
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for which we expect the relation Q < E < 1, is a measure 
for the "efficiency" of motion, i.e., Evq is the average 
speed at which the particles are able to move in their 
respective "target direction" e,. E i=s 1 corresponds to 
lanes, E — to a crystallized state. Representative sim- 
ulation results for the ensemble average (E) as a function 
of the noise intensity 9 are displayed in Figure E| 

We observed the following parameter dependencies: 
Crystallization is more pronounced for large r and large 
strip lengths L x , while small r and large strip widths L y 
are in favour of lane formation (but "freezing by heat- 
ing" still exists in the overdamped limit r — > 0). The 
number of particles required for crystallization does not 
depend on the length L x (if it is considerably larger than 
L y ), while it is roughly proportional to the width L y 
(for large enough L y ). Given a fixed aspect ratio L x /L y , 
for the system sizes that we could numerically handle 
there was no clear tendency whether the transition be- 
comes sharper or smoother with increasing system size 
N oc L x L y (see Fig. @). 

Another interesting quantity is the sum of the potential 
and kinetic energies associated with a given state: 



W= lim - 



dt 



Em o 
~9 V% 



-B 



(5) 



The paradox here is that the above mentioned crystal- 
lized state is usually more unstable than the fluid state, in 
the sense that the total energy (j^) of the system is higher 
in the crystallized state than in the "fluid" one (see inset 
of Fig. |^). Note that both, the "solid" (crystallized) state 
and the "fluid" state (i.e. lanes) are destabilized, if the 
friction term —vt/r is dropped, even in the case 9 = 
(see Ref. for a related inverse phenomenon). That 
is, without friction and due to the permanent driving, 
the undamped repulsive interactions eventually become 
destructive to any ordered state, which gives rise to a 
"gaseous" state. Therefore, we point out that the ener- 
getically less favourable crystallized state is maintained 
by the propulsion term m(w ej — Vi)/r in Eq. (^) which, 
by the way, is also relevant for lane formation. Note that 
the absolute value of this term becomes largest for Vi = 
(i.e. blocking), while it is small for the "fluid" state with 
Vi w v ei. 

We consider the transition to a stationary state with 
a higher total energy by increasing the noise intensity to 
be a signature of a novel class of behavior in certain non- 
equilibrium systems, which may have interesting appli- 
cations. However, here we have demonstrated "freezing 
by heating" only for limited sizes and a specific geome- 
try. Nevertheless, we point out that "freezing by heating" 
does not require walls, but occurs for periodic boundary 
conditions in y-direction as well. Also, we do not need 
periodic boundary conditions in x-direction. Any suffi- 



ciently long simulation area will produce both the orga- 
nized fluid state (i.e. lanes) and the crystallized state, if 
only the system is continuously entered by particles at 
the left-hand and right-hand boundaries. 

Why is freezing by heating new? Some glasses may 
crystallize when slowly heated. However, this is a 
well understood phenomenon: the amorphous state is 
metastable for those temperatures, and crystallization 
means an approach to the more stable state with smaller 
total energy. In general, one can distinguish three cases 
when fluctuations (i.e. temperature or external pertur- 
bations) are increased: (i) Total energy increases and 
order is destroyed (e.g., melting), (ii) Total energy de- 
creases, ordering takes place, and the system goes from 
a disordered metastable to an ordered stable state (e.g., 
in metallic glasses and some granular systems), (iii) To- 
tal energy increases and ordering takes place, while the 
system goes from a partially ordered stable to a highly 
ordered metastable state, which corresponds to the new 
situation presented here. 

In our case, crystallization is achieved by sponta- 
neously driving the system with the help of noise up- 
hill towards higher total energy. The system would like 
to maximize its efficiency H], but instead it ends up 
with minimal efficiency due to noise-induced crystalliza- 
tion. The role of "temperature" or noise here is to de- 
stroy the energetically more favourable fluid state, which 
inevitably leads to jamming and finally to crystall-likc 
lattices. The corresponding transition seems to be re- 
lated to the off-lattice nature of our model and is differ- 
ent from those reported for driven diffusive systems on a 
lattice PJlO). It should be noted that the transition we 
find is not sharp, which is a consequence partly of the 
mesoscopic nature of the phenomenon and partly of the 
disorder in the initial state. 

We would like to point out that "freezing by heating" 
is likely to be relevant to situations involving pedestrians 
under extreme conditions (panics). Imagine a very smoky 
situation, caused by a fire, in which people do not know 
which is the right way to escape. When panicking, people 
will just try to get ahead, with a reduced tendency to 
follow a certain direction. Thus, fluctuations will be very 
large, which can lead to fatal blockings. 

Our results demonstrate that in driven mesoscopic sys- 
tems phenomena qualitatively different from those occur- 
ring in thermodynamical systems can be observed. While 
most non-equilibrium transitions have analogies to equi- 
librium ones [l6| , the noise-induced ordering observed in 
the effect of "freezing by heating" is just opposite to the 
transitions occurring in equilibrium systems. This sug- 
gests that future studies along the lines of the present 
approach are likely to lead to further unexpected find- 
ings. 
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FIG. 1. Simulation of 20 particles moving from left to right 
(black) which interact with 20 particles moving from right to 
left (white) on a periodic strip of length L x = 20 and width 
L y = 5 at different noise intensities. The model parameters 
are m = 1, D = 1, v = 1, A = 0.2, B = 2, and r = 0.2. (a) 
Lanes of uniform directions of motion forming at small noise 
intensity (8 — 1). (b) Snapshot of an intermediate jammed 
state with a rough interface, which is about to form "chan- 
nels", (c) Final crystallized state resulting for large noise 
intensity (0 = 1000). 
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FIG. 2. Efficiency (E) of the system as a function of the 
noise intensity for various system sizes with a fixed aspect 
ratio of LxjL y = 4 : 1 and the parameter values displayed 
in Figure |l|. Shown are averages over 128 (for N = 160: 64) 
simulation runs with different random seeds. The inset shows 
the average potential, kinetic, and total energy in the result- 
ing system states. While the potential energy ( — ) increases 
with the noise amplitude due to jamming, the kinetic energy 

( ) is composed of a contribution proportional to and 

a decreasing contribution Nm(vo{E)) 2 jl. The total energy 
( — ) turns out to be an increasing function of 0. 



